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Renormalization-Group Evolution of the B-Meson 
Light-Cone Distribution Amplitude 



m 
O 
o 

(N 



> : 

<N . 
00 ■ 

O ■ 
m : 
o . 
m ■ 

o ; 

^ : 

a: 
i 

D 

^ ■ 



X 



Bjorn 0. Lange and Matthias Neubert 
Newman Laboratory for Elementary-Particle Physics, Cornell University, Ithaca, NY 14853, USA 

An integro-differential equation governing the evolution of the leading-order B-meson light-cone dis- 
tribution amplitude is derived. The anomalous dimension in this equation contains a logarithm of the 
renormalization scale, whose coefficient is identified with the cusp anomalous dimension of Wilson 
loops. The exact analytic solution of the evolution equation is obtained, from which the asymptotic 
behavior of the distribution amplitude is derived. These results can be used to resum Sudakov log- 
arithms entering the hard-scattering kernels in QCD factorization theorems for exclusive B decays. 
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1. Introduction. It has recently become apparent 
that the light-cone structure of the B meson is of great 
phenomenological interest. Using standard perturbative 
QCD methods for hard exclusive processes, B-meson 
light-cone distribution amplitudes (LCD As) were first in- 
troduced in a study of the asymptotic behavior of heavy- 
meson form factors at large momentum transfer [1]. Two 
such amplitudes, called <p±(w, /i), arise in the parameter- 
ization of -B-meson matrix elements of bilocal heavy-light 
current operators at leading-order in heavy-quark effec- 
tive theory (HQET). The interest in these quantities was 
revived when it was found that the amplitudes for many 
exclusive B-meson decays can be simplified significantly 
using QCD factorization theorems [2]. At leading power 
in Aqcd / irib these amplitudes can be expressed in terms 
of convolution integrals of perturbative hard-scattering 
kernels with the leading-order LCDAs. 

The £?-meson LCDAs appear in processes where light 
energetic particles are emitted into the final state, such 
as B — > 7T7T and B — > K*j, and at leading power arise 
from hard interactions of collinear partons with the soft 
spectator quark inside the B meson. In almost all appli- 
cations of QCD factorization theorems only the function 
(/>+(w,/i) contributes. Generically, the hard-spectator 
term in a QCD factorization formula is of the form 



du> 



(1) 



where the kernel T is calculable in perturbation theory. 
The factor X/uj ensures convergence of the integral for 
large uj. The characteristic scale for soft-collinear inter- 
actions is IEuj ~ m-fjAQCD, where E ~ m& is the energy 
of a collinear particle in the i?-meson rest frame. At fixed 
order in perturbation theory the kernel depends logarith- 
mically on this scale, and it is independent of u to lowest 
order in perturbation theory. In addition, the kernel can 
depend on hard scales of order nib- We will assume below 
that T is defined such that the convolution integral (1) 
is renormalization-group (RG) invariant. 

A thorough understanding of factorization requires 



controlling the scale dependence of the LCDA and the 
kernel under the convolution integral (1) using evolution 
equations. This is crucial for a clean separation of physics 
associated with different mass scales and for the system- 
atic resummation of large (Sudakov) logarithms, which 
enter the kernel at every order in perturbation theory. 
In this Letter we derive the RG equations for the LCDA 
</>+(w,/i) and for the kernel T(w,/i), present their exact 
analytic solutions, and extract model-independent results 
for the asymptotic behavior for u> — > and uj — > oo. 



2. Evolution equations. 
Fourier transform 
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The LCDA is given by the 



(2) 



of a function <p+(r, fi) defined in terms of a i?-meson ma- 
trix element in HQET. Denoting by h the effective heavy- 
quark field and by q s the soft spectator quark, and using 
a mass independent normalization of meson states, we 
write [1] 
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Here z || n with n 2 = is a null vector, v is the B- 
meson velocity, F represents an arbitrary Dirac matrix, 
and r = v ■ z — iO. The gauge string S n (z, 0) repre- 
sents a soft Wilson line connecting the points and z on 
a straight light-like segment. The quantity F(/i) is the 
HQET matrix element corresponding to the asymptotic 
value of the product f By/Trig in the heavy-quark limit. 
The distribution amplitude (f)^ (r, fj,) is normalized to 1 
at r = 0. The analytic properties of this function in the 
complex r plane imply that <^+(w, /i) = if oj < 0. 

We denote by 0+ (u>) the Fourier transform of the bilo- 
cal HQET operator in (3) and write the relation between 
bare and renormalized operators in the form 
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0™ n (uj,n) = J duj' Z+(u,u',n) 0^ c (uj'), 



(4) 



where Z + {uj,uj', /j,) = 5{u> — uj') at lowest order. Opera- 
tors with different momentum uj' can mix under renor- 
malization since they have the same quantum numbers. 
In the MS scheme the function Z + (uj, uj', /i) is defined so 
as to subtract the ultraviolet (UV) poles in the matrix 
elements of the bare operators. 

The B-meson matrix element of the renormalized op- 
erator T ? n (uj, n) is, up to a Dirac trace, given by the 
product F(n) </>+(w, /x). It follows that the LCDA obeys 
the evolution equation 



J^^+i^^) = -J cL;'7 + (u;,u/,/i)<^?(u/,/x) (5) 

with the anomalous dimension (unless otherwise indi- 
cated, a s = a s (n)) 



7 + (w,u/,/j) 



/ 



_ dZ + (u>, uj, /i) 



-1f(<u s )5{uj - J) . 



Z+ (d>,w',/i) 



(6) 



Here 7f is the universal anomalous dimension of lo- 
cal heavy-light currents in HQET, which determines the 
scale dependence of F(fi). 

At one-loop order the result for Z + (uj,uj', /j) is ob- 
tained by evaluating the UV poles of the one-loop HQET 
diagrams shown in Figure 1. The loop integrals con- 
tain 5-functions whose arguments depend on the value 
of the "plus component" k + = n ■ k of the loop mo- 
mentum. The integration over the orthogonal "minus 
component" is most easily performed using the theorem 
of residues, and the remaining integration over trans- 
verse components can be evaluated using standard tech- 
niques of dimensional regularization (setting d± = 2 — 2e 
for the dimension of the transverse space). In the first 
two diagrams the transverse loop integrations give rise 
to UV poles, whereas the third diagram is UV finite. 
In units of CfOLs/^-k, the corresponding contributions to 
Z + (uj,uj', n) are given by (k = — k + ) 



dk 
~k 



D, = — 



-2e 



5{k + uj' - uj)- 6(uj - uj') 



dk uj' 
k u 



(7) 



5{k + uj-uj') -5(uj-uj') 



These expressions must be understood in the sense of 
distributions, and only pole terms in 1/e are kept after 
integration over uj' in (4) . Note that the k- integral in the 
first diagram in Figure 1 gives rise to a 1/e 2 pole times 
S(uj — uj'). Accounting for wave-function renormalization 
of the external quark fields we obtain 



Z ( + 1) (uj,uj'^) = 



1 
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^ S(uj - uj') 



uj 6(uj'-uj) 6(uj-oj') 
uj' uj' — UJ uj — uj' 



(8) 
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FIG. 1. One- loop diagrams contributing to the calculation 
of the anomalous dimension. The crossed circle denotes an 
insertion of the operator 0+ arc (a/). The double lines represent 
heavy-quark fields in HQET. 



Here and below a superscript "(1)" is used to indicate 
one-loop coefficients in units of CfolsI^-k. The plus dis- 
tribution is defined such that, when Z + is integrated with 
a function f(uj'), one must replace f(uj') — > f(uj') — ,f(uj) 
under the integral. The non-diagonal terms in (8) agree 
with those found in [1]; however, the double pole and log- 
arithm in the first term were missed in that paper. For 
the anomalous dimension in (6) we now obtain 



7 + (w,a/, / u) = r cuS p(a s ) In- +7(0;^) 5{uj - uj') 
+ ujT(uj,uj', a s ) , 



(9) 



where / duj' T(uj, uj' , a s ) — by definition. Taking into 

account that 7^ = —3 is the one-loop coefficient of the 
anomalous dimension of heavy-light currents, it follows 
from (8) that T^lp = 4, 7W = -2, and 



r (1) (^') = -r«I» P 



6(uj' -uj) 9(uj - uj') 



uj'(uj' — uj) uj{uj — UJ 1 ) 



(10) 



The anomalous dimension in (9) is of the Sudakov type. 
It contains a logarithmic dependence on the renormal- 
ization scale in addition to the dependence through the 
running coupling a s {p). This feature distinguishes the 
kernel for the B- meson LCDA from the familiar Brodsky 
Lepage kernel for the LCDA of a light pscudoscalar me- 
son [3]. The extra logarithm has its origin in the renor- 
malization properties of Wilson lines with light-like seg- 
ments [4]. Using the fact that a heavy quark in HQET 
can be described by a Wilson line h(0) = S v (0, —00) ho, 
where ho is a sterile field without QCD interactions, it 
follows that the interacting fields in the operator in (3) 
can be written as q s (z) S n (z, 0) S v (0, — 00). This corre- 
sponds to a gauge string extending from minus infin- 
ity to along the v direction, another string extending 
from to z along the light-like n direction, and a light 
quark field located at the end of the string at position 
z. This Wilson line has a cusp singularity at the origin, 
which gives rise to a local, single-logarithmic term in the 
anomalous dimension. From (7) it is seen that the log- 
arithm arises indeed from the gluon exchange between 
the two strings S n and S v (first diagram in Figure 1). 
Such cusp singularities are universal and of geometric 
origin. The corresponding anomalous dimension r cusp is 
process independent and has been computed at two-loop 
order in [4]. The important observation following from 
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this discussion is that also in higher orders there is only 
a single logarithm present in the anomalous dimension 
7+(w, u)',/j,). With this knowledge, the evolution equa- 
tion can be solved in RG-improved perturbation theory. 

From the scale independence of the convolution inte- 
gral (1) it follows that the hard-scattering kernel obeys 
the evolution equation 



the exact expressions in (15) one can derive approximate 
results for the functions g and U at a given order in RG- 
improved perturbation theory. The explicit forms arising 
at next-to-leading order can be found in [5] . 

We now assume that the function </>+(w, no) is given at 
some low scale [1q ~ few x Aqcd arid define its Fourier 
transform with respect to ln(w/A*o) through 



r°° u) i r°° f 1 1 \ u 

T{u,ti)= dw , -7+(w',w,/i)T(w', A1 ). (11) 0?(w,/io) = ^-/ dt<p (t){ — ) > 



(16) 



The explicit result for the anomalous dimension given 
above shows that ^y^ + (w',u), fi) = , fi) (at least 

to one- loop order). Hence, apart from an overall sign the 
LCDA and the kernel obey the same RG equation. 

3. Analytic solutions. The general solution of the 
evolution equations (5) and (11) can be obtained using 
the fact that on dimensional grounds 



/ 



duj' cj T(w, w', a s ) (uj') a = u a 5{a, a s ) , 



(12) 



where the function "5 only depends on the exponent a and 
the coupling constant, and 3^(0, a s ) — by definition. 
The integral on the left-hand side is convergent as long 
as —1 < Re a < 1. The corresponding integral with 
T(w, w', a s ) replaced by F(o/, u>, a s ), which is relevant to 
the evolution equation (11) for the kernel, is given by 
3 r (— a, a s ). This follows from the fact that u', a s ) — 
uj~ 2 f(ui'/u)) on dimensional grounds. At one-loop order 
we find from (10) 

3™ (a) = r« p [V(l + a) + ^(1 - a) + 2 lE ] , (13) 

where ip(z) is the logarithmic derivative of the Euler T- 
function. 

Relation (12) implies that the ansatz 



with 

g(fj,,Ho) = 



-) 

Mo/ 



expU(fj,,fj, ,v) ( 14 ) 



da rcu5!> 



(a) 



(15) 



U(fi,no,r)) 



da 



g{li, fj, a ) + 7(a) 



provides a solution to the evolution equation (5) with ini- 
tial condition f(w, fio, Mo, v) — ( w // u o)' ! at some scale /xo- 
Here /i a is defined such that a s (fj, a ) = a, the /3-function 
is P(a s ) — da s /d\n.n, and r] can be an arbitrary com- 
plex parameter. Note that g([i, /j®) > if (i > Mo- Given 



where <po{0) = 1/\b is determined in terms of the first 
inverse moment of the LCDA at the scale Mo [2] . It follows 
that the result for the LCDA at a different scale \i is given 

by 

l f°° 

A«) = ^ J dtipo(t) f(u,(j,,fio,it) . (17) 

It is instructive to work out the t and w-dependence of the 
integrand at leading order in RG-improved perturbation 
theory. Using the one-loop result (13) we obtain 



/(<*>, M> Mo,**) oc 



it+g 



T(l - it - g)T{\ + it) 
T(l + it + g)T(l-it) 



, (18) 



where g = (2Cf//?o) ln[a s (iJ.o)/ a s(fj)] is the leading- 
order contribution to the function g(/j,,fi ), and the for- 
mula can only be trusted if g < 1 (which is justified for 
all reasonable parameter values). The expression above 
is an analytic function in the complex i-plane with singu- 
larities along the imaginary axis. For \i > fi the nearest 
singularities are located at t = — g) and t = i. 

The locations of the nearest singularities in the product 
ifo(t) /(w, m, Mo, it) determine the asymptotic behavior of 
the LCDA for uo — > (lower half-plane) and uo — > oo (up- 
per half-plane). If we assume that the function <fi+(uj, Mo) 
at the low scale vanishes like uu s for small u> and falls off 
like for large uu (exponential fall-off would correspond 
to £ — > oo), then the nearest poles in ifo{t) are located 
at t = —i5 and t — i£, respectively. It follows that after 
evolution effects 



(j>+{u,n) 



min(l, S-\-g) . 



for w — > 0, 

w - min(l,€)+9 . forw^OO. 



(19) 



Irrespective of the initial behavior of the LCDA, evolu- 
tion effects drive it toward a linear growth at the origin 
and generate a radiative tail that falls off slower than 1/uj 
even if the initial function has an arbitrarily rapid fall- 
off. This implies, in particular, that the normalization 
integral of (f>?_(u), m) is UV divergent for large values of w. 
This fact, which was already noted in [1], is not an ob- 
stacle to our analysis. Convolution integrals appearing in 
QCD factorization theorems are always of the form (1), in 
which the extra 1/u) factor suppresses the contributions 
from large w and renders the integral finite. 

The solution of the evolution equation (11) for the 
hard-scattering kernel proceeds in a similar way, ex- 
cept that in the expression for the function / in (14) 
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one must replace jio — > /Uj, g(jj,,fj,o) — > — giPiHi), an d 
[/(/x, /Lto, r/) — ► —U(fi, Hi, —if). As mentioned in the Intro- 
duction, at fixed order in perturbation theory and at an 
intermediate scale Lii ~ V TO &A, the kernel depends on o> 
only through logarithms of the type ln(2£ , w//z?), which 
are of order 1 and so do not need to be resummed. We 
can therefore write T(u),/j,i) = 7[\n(2Euj / Lif ), . . .], where 
the dots represent other arguments independent of lo. 
(Some of these arguments may contain large logarithms, 
which must be resummed separately.) It follows that by 
considering derivatives of the solution with respect to i] 
(evaluated at r\ = 0) we can satisfy arbitrary initial condi- 
tions at the scale Lii- We can then solve (11) and compute 
the hard-scattering kernel at a scale \x < iii. The exact 
solution is given by 



T(w, M ) = T[V r) ,...] 



(2Ew\ n 



\ A ) VM 
x exp[-U(n,iJ,i,-r])} 



9(w,fi) 



?7=0 



(20) 



where the notation T[V,,, . . .] means that one must re- 
place each logarithm of the ratio 2Eui/Lij in the ini- 
tial condition by a derivative with respect to the aux- 
iliary parameter n. It follows that the kernel scales like 
T <~ modulo logarithms. 

From the explicit solutions given above it can be seen 
that the convolution integral in (1) is indeed indepen- 
dent of the renormalization scale. For this, it is nec- 
essary to move the integration contour in (17) into the 
upper complex t plane by an amount 77 + g(iii, liq). Note, 
however, that the product T(lo, /j,) 4>+(uj, ii) for fixed lo 
is not scale independent. We also stress that evolution 
effects mix different (fractional) moments of the LCDA. 
For instance, the first inverse moment of the LCDA at 
a scale fi, which is sometimes called 1/A_b(/x), is related 
to a fractional inverse moment of order 1 — g{fi, y!) at a 
different scale it! . As a result, the scale dependence of the 
parameter A b (p) is not calculable in perturbation theory. 
Controlling it would require knowledge of the functional 
form of the LCDA. 

To illustrate our results we consider a scenario where 
the hard-scattering kernel at a high scale li% with 
OL s {\ii) — 0.3 takes the simple form T(ui,iii) = 1, 
and where the LCDA at a low hadronic scale 
with a s (fio) = 1 assumes the form <P+(lo,lio) = 
(lo/X b ) e -^/ x B suggested in [1], for which 



(21) 



Figure 2 shows the results for the LCDA and the ker- 
nel at three different values of the renormalization scale. 
For simplicity, the anomalous dimensions and /3-function 
are evaluated at one-loop order. Whereas the kernel ex- 
hibits a smooth power-like behavior, the most character- 
istic feature of the evolution of the LCDA is the devel- 
opment of a radiative tail for large values of ui. In this 




1 1.5 2 2.5 

lo [GeV] 

FIG. 2. Results for the LCDA and the kernel (inset) 
for different values of the renormalization scale such that 
a 3 = 1 (solid), 0.5 (dashed), 0.3 (dashed-dotted) . We use 
A_g = 0.3 GeV at fi — no. 



example the value of the convolution integral (1) is about 
25% smaller than the value 1/A_b which one would obtain 
without evolution effects. 

4- Conclusions. We have derived evolution equations 
for the leading-order B-meson LCDA and for the hard- 
scattering kernels entering the spectator term in QCD 
factorization theorems for exclusive B decays into light 
particles. Simple scaling relations are obtained for the 
asymptotic behavior of these quantities. Exact analytic 
solutions to the evolution equations are given in terms 
of integrals over anomalous dimension functions. This 
accomplishes the resummation of large Sudakov double 
logarithms to all orders in perturbation theory. 
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